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Wolf-Rayet (WR) stars, as they are advanced stages of the life of massive stars, provide a good 
test for various physical processes involved in the modelling of massive stars, such as rotation 
and mass loss. In this paper, we show the outputs of the latest grids of single massive stars 
computed with the Geneva stellar evolution code, and compare them with some observations. 
We present a short discussion on the shortcomings of single stars models and we also briefly 
discuss the impact of binarity on the WR populations. 


1 Introduction 

Wolf-Rayet (WR) stars are known to be the very 
advanced stages of the life of massive stars (Conti 
1976). These stars exhibit at their surface signs of 
hydrogen burning products through CNO cycle (for 
the WN stars) or even of helium burning products 
(Crowther 2007). To explain the existence of such 
objects, elements that are produced in the centre of a 
massive stars have to appear at the surface. This can 
be achieve by two different ways: a) internal mixing 
inside that star, allowing for transporting chemical 
elements from the place where they are produced to 
the surface, b) mass loss, that progressively remove 
the external layer of the star, and can ultimately un¬ 
cover deep regions of the star, where nuclear burning 
previously occurred. 

Internal mixing can be induced by the convective 
movements of matter in the convective regions, or 
by any other kind of mixing process in the radiative 
regions, such as rotational mixing (e.g. Zahn 1992). 
Mass loss can be due either by stellar winds, or by 
mass transfer in close binaries. 

In this paper, we present preliminary results show¬ 
ing what are the expected WR stars population in 
the single massive star scenario, at different metal- 
licities. The assumptions made in our computations 
and the chosen prescriptions are shortly reminded in 
section 2. Then, the various mass limits and num¬ 
ber ratio of WR subtypes are presented in section 3. 
Finally, we briefly discuss the interesting challenges 
that need to be resolved to make progresses in that 
field. 

2 Ingredients of the stellar 
models 

The results presented in this paper were obtained 
by analysing the models from Ekstrom et al. (2012), 
Georgy et al. (2013), and Eggenberger et al. (in 
prep.). The computations were performed using the 
Geneva stellar evolution code (Eggenberger et al. 
2008). Here is a brief summary of the assumed phys¬ 
ical ingredients and prescriptions we used. 

Convective zones are determined with the 


Schwarzschild criterion. For hydrogen and helium 
burning cores, a penetrative overshoot is assumed, 
which extends over a fraction dover of the pressure 
scale height at the formal edge of the convective core. 
Inside the convective zones, the mixing is assumed 
to be instantaneous. In the inner regions, convec¬ 
tion is assumed to be adiabatic, while in the ex¬ 
ternal envelope, the thermal gradient is computed 
in the framework of the mixing-length theory (MLT 
Bohm-Vitense 1958). The free parameter dover is 
calibrated at solar metallicity, by reproducing the 
width of the main sequence (MS) at around 1.7 Mq 
(Ekstrom et al. 2012). The value of the a parameter 
appearing in the MLT is calibrated on a solar model. 

Rotation is included in the framework of the “shel- 
lular” rotation approximation (Zahn 1992; Meynet 
& Maeder 1997; Maeder & Zahn 1998). The trans¬ 
port of angular momentum is computed with both 
the advective and diffusive term (Zahn 1992), while 
the transport of chemical elements is a purely dif¬ 
fusive process (Chaboyer & Zahn 1992). In this 
context, two diffusion coefficients are required. In 
our works, we use Di, from Zahn (1992) and Ughear 
from Maeder (1997). This choice best reproduces 
the presence of Gepheid loops at solar metallicity. 
The free parameter appearing in the formulation of 
7?shear is Calibrated in order to reproduce the ob¬ 
served enrichment of stars in the range of 10 to 
20 Mq at solar metallicity (Ekstrom et al. 2012). 

The mass-loss prescriptions and metallicity depen¬ 
dence are described in Ekstrom et al. (2012) and 
Georgy et al. (2013). An enhanced mass-loss rate 
(by a factor of 3) is used for red supergiant stars 
(RSG) above 20 Mq (see discussion in Georgy 2012; 
Georgy et al. 2012). 

In this work, we use the following criteria to clas¬ 
sify the WR stars (Georgy et al. 2012). We con¬ 
sider that a star is a WR stars when its log(reff) 
is above 4., and its surface mass fraction is smaller 
than 0.3. A WR star with more than 10“^ as sur¬ 
face H mass fraction is a WNL star. If the surface 
hydrogen abundances drops below 10“®, we have a 
WNE star if the surface nitrogen abundance is big¬ 
ger than the carbon one, and a WC star otherwise. 
These are obviously rough criteria, but they have 
the advantage that they can be computed with stan¬ 
dard outputs from stellar evolution models. A more 


1 



C. Georgy et al. 





Fig. 1: HRD of rotating massive star models: Z — 0.002 (SMC, kft panel), Z = 0.006 (LMC, middle panel), 
and Z = 0.014 (solar, right pannel). In each panel, models from 15 Mq (bottom track) to 120 Mq (top track) are 
shown. The different WR stages are indicated in colour: WNL (green), WNE (purple) and WC (blue). The initial 
rotation velocity is 40% of the critical one. Models are taken from Georgy et al. (2012) for Z = 0.002, Eggenberger 
et al. (in prep.) for Z = 0.006 and Ekstrom et al. (2012) for Z = 0.014. 


detailed classification would require the modelling 
of the emergent spectrum (taking into account the 
winds), and this is so far not possible during the 
computation runtime. Such spectrum modelling in 
a post-processing approach has shown that a classi¬ 
fication based on the aspect of the spectrum leads 
to significant differences (Groh et al. 2014). 

3 Wolf-Rayet stars from single 
star models 

Figure 1 shows the Hertzsprung-Russell diagram 
(HRD) for our latest models of massive stars, for 
15, 20, 25, 32, 40, 60, 85, and 120 Mq, at three 
metallicities: Z = 0.002 (left panel), corresponding 
to the metallicity of the Small Magellanic Cloud, 
Z = 0.006 {middle panel), corresponding to the 
Large Magellanic Cloud one, and finally solar metal¬ 
licity Z = 0.014 {right panel). The successive WR 
stages are highlighted with different colours: the 
WNL phase is shown in green, the WNE one in pur¬ 
ple, and the WC one in blue. 

A direct consequence of the weakness of the stel¬ 
lar winds at lower metallicity (Vink et al. 2001; El- 
dridge & Vink 2006) is that it becomes more and 
more difficult to produce WR stars when metallic¬ 
ity is decreasing. Moreover, at solar metallicity, our 
models spend a longer time in the RSC phase, allow¬ 
ing them to lose a lot of mass, and helping them to 
cross the HRD towards the WR phase at lower mass 
(around 20 Mq, see also Vanbeveren et al. 1998). 

The minimal masses to enter into a given WR 
phase are given in Table 1. At solar and LMC metal¬ 
licities, all WR types are accessible to single star 
models. Note however that the WC subtype oc¬ 
curs only for very massive stars at the LMC metal¬ 
licity, while at solar Z, models more massive that 


about 27 Mq are able to reach this stage. More¬ 
over, these stars are found at luminosities higher 
than log(L) ~ 5.3 at Zq and log(L) ~ 6.2 at Zlmc- 
At the metallicity of the SMC, the radiative stellar 
winds are too weak to produce stars more evolved 
than the WNL phase, and no stars without hydro¬ 
gen at the surface are expected from our models. 
The minimal mass to produce a WR star progres¬ 
sively shift from about 20 Mq at Zq to more than 
50 Mq at the SMC metallicity. Above these masses, 
the endpoint of the evolution of our models are thus 
WR stars from our simple classification scheme (see 
however Croh et al. 2013b,c). It is not yet clear what 
would be the final fate of this kind of stars in terms 
of supernova explosion (type Ibc, failed supernova, 
direct collapse? See Heger et al. 2003; Dessart et al. 
2011; Croh et al. 2013a; Bersten et al. 2014). 



0-star 

WNL 

WNE 

WC 

Solar 

15.8 Mq 

20.0 Mq 

25.3 Mq 

27.0 Mq 

LMC 

14.2 Mq 

32.1 Mq 

60.8 Mq 

63.1 Mq 

SMC 

12.6 Mq 

53.5 Mq 

- 

- 


Tab. 1: Minimal mass to enter into a given phase. Re¬ 
sults at Zq are taken from Georgy et al. (2012). 

The expected ratio of WR stars to 0-type stars, 
as well as subtype ratios, are shown in Table 2, in 
the constant star formation context. These num¬ 
bers confirm that the number of WR stars with re¬ 
spect to 0-type stars is decreasing with decreasing 
metallicity. The fraction of WC star is also ex¬ 
pected to decrease at low metallicity, as well as the 
WC/WN ratio. Comparison with the observation 
at solar metallicity was presented in Ceorgy et al. 
(2012). Our results showed that in order to repro¬ 
duce the observed ratio, WR stars should originate 
from a stellar population containing about 50% of 
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WR/O-stars 

WNL/WR 

WNE/WR 

WN/WR 

WG/WR 

WC/WN 

Solar 

0.066 

0.687 

0.022 

0.709 

0.291 

0.409 

LMG 

0.016 

0.887 

0.005 

0.892 

0.108 

0.121 

SMG 

0.006 

1.000 

0.000 

0.000 

0.000 

0.000 


Tab. 2: Various type ratios. Values at Zq are taken from Georgy et al. (2012). 


binaries, in good agreement with results from more 
elaborated synthetic binary population codes (Van- 
beveren et al. 1998; Eldridge et al. 2008). Note also 
that our WC/WN ratio reproduce reasonably well 
the observed trend at different metallicities (Neu- 
gent et al. 2012). 

4 Discussion 

Since a few years, it became clear that massive star 
populations contain a significant fraction of binary 
stars (e.g. Sana et al. 2012). It is thus important 
to know what are the physical mechanisms respon¬ 
sible for the appearance of the WR phenomenon (or 
subtypes), and what is linked to binary evolution or 
not. Recent observations of Galactic WC stars show 
a significant number of such kind of stars at surpris¬ 
ingly low luminosities (log(L) ~ 5.2, Sander et al. 
2012). These stars are hard to form through sin¬ 
gle star channel, even with a strong mass-loss rate 
during the RSG phase (Georgy 2012; Meynet et al. 
2015), except if we assume that the strong mass loss 
continues after the RSG phase (Vanbeveren et al. 
1998). On the other hand, they are routinely pro¬ 
duced through the binary channel (Eldridge et al. 
2008, 2013). 

The relatively high number of high luminosity 
WNL stars with a large fraction of hydrogen on their 
surface (e.g. Hainich et al. 2014) on the other hand 
points to the need of accounting for a proper treat¬ 
ment of the internal mixing of the star, particularly 
inside the radiative zones (Georgy et al. 2012). This 
illustrates the need of a correct treatment of the 
physics of single massive star, that definitely inter¬ 
venes as well in the modelling of binary evolution. 
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